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Abstract 

Measurement of the T{h — > 77)Br(/i bb) at the Photon Colhder at TESLA is studied 
for the Standard Model Higgs boson with mass of 120 to 160 GeV. The NLO estimation of 
background, analysis of overlaying events, realistic 6-tagging and corrections for escaping 
neutrinos were performed. We find that for Mh = 120-160 GeV the T(h — > 77)Br(/i — > bb) 
can be measured with a statistical accuracy of 2-7% after one year of the Photon Collider 
running. 
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1 Introduction 



A search of the Higgs boson is among the most important tasks for the present and future 
coUiders. Once the Higgs boson is discovered, it will be crucial to determine its properties with 
a high accuracy. A photon-collider option of the TESLA collider [1] offers a unique possibility 
to produce the Higgs boson as an s-channel resonance. The neutral Higgs boson couples to the 
photons through a loop with the massive charged particles. This loop-induced hjj coupling is 
sensitive to contributions of new particles which appear in various extensions of the SM. 

The SM Higgs boson with a mass below ~ 140 GeV is expected to decay predominantly 
into the bb final state. Here we consider the process ^ h ^ bb for a Higgs-boson mass 
Mh = 120, 130, 140, 150 and 160 GeV at the Photon CoUider at TESLA. Both the signal and 
background events are generated according to a realistic photon-photon luminosity spectrum 
[2], parametrized by a CompAZ model [3]. For the first time in such study overlaying events 
77 —>■ hadrons are taken into account, for which we use photon-photon luminosity spectra from 
a full simulation [2] . Our analysis incorporates a simulation of the detector response according 
to the program SIMDET [4] and - the next new element - a reahstic 6-tagging [5]. This analysis 
supersedes our earlier analyses presented in [6, 7, 8]. 



2 Photon— photon luminosity spectra 

The Compton back-scattering of a laser light off high-energy electron beams is considered as a 
source of high energy, highly polarized photon beams [9] . According to the current design [1] , 
the energy of the laser photons is assumed to be fixed for all considered electron-beam energies; 
laser photons are assumed to have circular polarization Pc — 100%, electrons longitudinal 
polarization is Pg = 85%. We use the luminosity spectrum peaked at high energy and assume 
that the energy of primary electrons is adjusted in order to enhance the signal at a particular 
mass. 

In a generation of the processes 77 — > hadrons one has to take into account also the low 
energy events, since they contribute to overlaying events [1]. To simulate them we use the 
realistic 77 luminosity-spectra for the photon collider at TESLA [2], with the non- linear cor- 
rections and higher order QED processes. For generation of the processes 77 — > /i — > 66 and 
77 — > bb{g), cc{g) we use the CompAZ parametrization [3] of the spectrum [2]. 

The results presented in this paper were obtained for an integrated luminosity corresponding 
to one year of the Photon Collider running, as given by [2]. For example, for y/s^ — 210.5 GeV, 
which is optimal for Mh = 120 GeV, the total photon-photon luminosity per year is L^^ = 410 
fb~^ (84 fb~^ for Wyy > 80 GeV). The total photon-photon luminosity increases to about 490 
fb-^ for ^ = 260 GeV (used for Mh = 160 GeV). 
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Transverse energy per bunch crossing 




Figure 1: Angular distributions of transverse energy, Et, for 77 hadrons events per bunch 
crossing. Various components and their sums are indicated. Generation was done for -y/s^ — 
210.5 GeV [8]. 

3 Details of a simulation 

and the first results for Mh =120 GeV 



We calculated the total width and branching ratios of the SM Higgs boson, using the program 
HDECAY [10], where higher order QCD corrections are included. A generation of events 
was done with the PYTHIA 6.214 program [11]. A parton shower algorithm, implemented in 
PYTHIA, was used to generate the final-state particles. 

The background events due to processes 77 — > bb{g), cc{g) were generated using the program 
written by G. Jikia [12], where a complete NLO QCD calculation for the production of massive 
quarks is performed within the massive-quark scheme. The program includes exact one-loop 
QCD corrections to the lowest order (LO) process 77 —>■ bb, cc [13], and in addition the non- 
Sudakov form factor in the double-logarithmic approximation, calculated up to four loops [14]. 

For an estimation of systematic uncertainties in 6-tagging simulation we also generated the 
background events, using the LO QED cross section for the processes 77 —>■ bb and 77 cc, 
and including parton shower, as implemented in PYTHIA. 

The fragmentation into hadrons for all processes was performed using the PYTHIA program. 
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Because of the large cross section, about one 77 — ^ hadrons event^ is expected per bunch 
crossing at TESLA Photon CoUider (for ^/s^ = 210-260 GeV, with nominal luminosity). We 
generate these events according to PYTHIA 6.214, including direct and hadron-hke photon 
contributions. We use the full simulation of photon-photon spectra [2], rescaled to the chosen 
beam energy. For each considered collision energy, ^/s^, average number of 77 hadrons 
events per bunch crossing is calculated. Next, for each signal '-ff ^ h ^ bb or background 
77 —>■ bb{g), cc{g) event 77 — >• hadrons eyents are overlaid (added to the event record) according 
to the Poisson distribution. 

Processes contributing to the overlaying events have forward-peaked distributions, as is 
shown in Fig. 1 for ^/s^ = 210.5 GeV. Therefore, to minimize an influence of these events 
on the signal measurement we ignore tracks and clusters with |cos(^i)| > cos(^min) — 0.9 
(dmin = 450 mrad; where the angle between the beam axis and a track/cluster, 6i, is measured 
in the laboratory frame). Hadron-level studies performed with PYTHIA show that after this 
cut overlaying events contribute to the measured invariant mass below 5% in 90% of signal 
events. Below we use the Omin cut only when overlaying events are included in the analysis. 

The fast simulation program SIMDET version 4.01 [4] was used to model a TESLA detector 
performance. 

Jets were reconstructed using the Durham algorithm, with ycut = 0.02; the distance measure 
was defined as yij = 2mm{E^, Ej){l — cos9ij)/E^-g, where E^is is the total energy measured in 
the detector. 

The following cuts were used to select the h ^ bb events: 

1. since the Higgs bosons are expected to be produced almost at rest, we require that the 
ratio of the total longitudinal momentum of all observed particles to the total visible 
energy is \P^\/Eyis < 0.15, 

2. we select two- and three-jet events, Nj^ts = 2, 3, so that events with one additional jet 
due to a hard-gluon emission are also accepted, 

3. for each jet we require | cos^i| < 0.75, i = 1, Njets- 

We use "ZVTOP-B-Hadron- Tagger" package for the TESLA collider [5] for realistic b- 
tagging simulation. The package is based on the neural-network algorithm trained on the 
Z decays. For each jet it returns a "6-tag" value - the number between and 1 corre- 
sponding to "6-jet" likelihood. At y^s^ = 210.5 GeV, the signal to the background ratio, 
A^(77 —> bb{g))/N{'-f'-f cc{g)), for selected 2-jet events, after additional cut E^s > 85 GeV, 
was investigated as a function of two 6-tag values, with and without overlaying events. 

By accepting 77 — > bb{g), cc{g) events above a given signal-to-background ratio, the bb- 
tagging efficiency Ebb was studied as a function of cc-mistagging probability Ecc, as shown in 
Fig. 2. For 3-jet events three possible pairs of jets were considered and the event was accepted 

^We consider only photon-photon events with W^^ > 4 GeV. 
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if at least one pair gives the signal-to-background ratio above the cut. It was found that the cut 
corresponding to the efficiencies Ebb = 80% and Ecc = 2.2% is optimal for the T{h — > 77)Br(/i 
bb) measurement. For other considered electron-beam energies similar efficiencies were obtained 
(not shown). In the earlier analyses [12, 6] a fixed efficiency for the 66-tagging, Ebb = 70%, and 
a fixed probability for a mistagging of the cc events, Ecc = 3.5%, were assumed (indicated in 
Fig. 2). 




Figure 2: The 66-tagging efficiency of 77 bb{g) events. Ebb, versus cc-mistagging probability 
of 77 — > cc{g) events, Ecc, for ^/s^ = 210.5 GeV with the additional cut Ems > 85 GeV, without 
and with overlaying events. Optimal Ebb (and Ecc) from these simulations (square and dot) and 
the earlier estimate (star) are indicated. 

In order to estimate possible influence of neglecting soft-gluon emissions in simulation of a 
NLO background, for which parton shower algorithm is not used, the efficiency of 6&-tagging 
was compared with an efficiency for the LO background simulation with PYTHIA (including 
parton shower) [7]. The overlaying events were not included in this investigation. Results of 
this analysis do not indicate any significant influence on the 6-tagging. We obtain efficiencies 
of Ebb = 82.2% (82.4%) and Ecc = 2.2% (2.3%), for background events generated according to 
NLO cross-sections without parton shower (LO cross-sections with parton shower). 

In Fig. 3 we show influence of a 6-tagging on reconstructed invariant-mass distribution, Wrec, 
for the signal events 77 /i — bb. These results were obtained without overlaying events 
77 hadrons. The low mass tail is due to the presence of events with escaping neutrinos 
(see [6] for more details). Contribution of these events can be suppressed by an additional cut 



5 



e e beams with Vs„„ =210.5 GeV 



> 2250 

o 

^ 2000 

(» 

a 1750 
Qi 

Qi 

1500 

Qi 

I 1250 

s 

1000 



750 



500 



250 



Before b-tagging: 

no Pj/Ej cut 

• • • • p^^cut 

After b-tagging: 
- - - - no p.p/E.p cut 
p.p/E,p cut 



= 120 GeV 
Total = 410 fb"^ 




140 160 

W_ (GeV) 



Figure 3: Reconstructed invariant mass, Wrec, distributions for selected —>■ h bb events, 
for Mh = 120 GeV. Distributions obtained before and after 6-tagging, without and with an 
additional Pt/Et < 0.04 cut are compared. 



Prp/Erp < 0.04, where Pt and are the absolute values of the total transverse momentum 
of an event, Pt, and the total transverse energy, respectively. We see that 6-tagging does not 
influence significantly the shape of the distributions. 

In Fig. 4 we compare the invariant- mass distributions before and after taking into account 
the overlaying events. A mass resolution, derived from the Gaussian fit in the region from fi — a 
to fi + 2a, is 3.5 and 6.1 GeV, respectively. Despite a quite high 9min cut, the overlaying events 
result in a bigger tail above Wrec = 120 GeV. A small drop in a selection efficiency, resulting in 
the reduced number of events (from about 8520 to 7740 events), is observed. This is because 
the energy deposits from the 77 — > hadrons processes, remaining after the Omin cut, "shift" 
jets nearer to the beam axis and the event can be rejected by the jet-angle cut. Moreover, the 
additional deposits and ^mm-cut deform jets, what slightly reduces the selection efficiency. To 
study this issue in details we plan to simulate in the future the signal events with various Omin 
values. 

After applying the selection cuts, 6-tagging and rejecting low-angle deposits we obtain the 
distributions of the reconstructed 77 invariant mass, Wrec-, shown in Fig. 5. The signal and 
NLO background contributions, hb{g) and cc{g), are shown separately. 

Assuming that the signal for Higgs-boson production will be extracted by counting the 
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Figure 4: Reconstructed invariant-mass, Wrec, distributions for selected •yy h ^ bb events 
after 6-tagging, for = 120 GeV, obtained without and with overlaying events (OE). 

number of bb events in the mass window around the peak. Nobs, and subtracting the expected 
background events, Nbkgd, we can calculate the expected relative statistical error for the cross 
section (7(77 ^ h ^ bb) (or the partial width multiplied by the branching ratio r{h 
77)Br(/i — i> bb)) in the following way: 



A(t(77 ^h^bb) A [r(/i ^ 77)Br(/i ^ bb)\ ^ 



C7(77 ^ h^bb) r{h 77)Br(/i bb) 



Nobs — Nbkgd 



The accuracy expected for the Higgs-boson mass of 120 GeV, from the reconstructed invariant- 
mass distribution in the selected mass range between 102.5 and 142.5 GeV (see Fig. 5), is equal 
to 2.0%. It is in agreement with the results of a previous analysis [6]. Note however, that in 
the present analysis a loss of signal efficiency due to the overlaying events is compensated by a 
more effective 6-tagging. 



4 Final results for masses 120—160 GeV 

As in [6], to correct for escaping neutrinos we use the corrected invariant mass, a variable 
defined as: 



Wcorr = V^rec + 2Pt(^«s + Pt) (1) 
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Figure 5: Reconstructed invariant-mass, Wrec, distributions for selected bb events. Contribu- 
tions of the signal, for = 120 GeV, and of the heavy-quark background, calculated in the 
NLO QCD, are shown separately. Arrows indicate the mass window, 102.5 to 142.5 GeV, opti- 
mized for the measurement of the r{h — > 77)Br(/i —>■ bb), which leads to the statistical precision 
of 2.0%. 

In Fig. 6 the distributions of Wcorr for the selected signal events, without and with overlaying 
events, are presented. The tail of events with invariant masses below ~ 110 GeV is much smaller 
than for the VFrec-distributions (compare with Fig. 4). The mass resolutions, derived from the 
Gaussian fits to the VFcorr-distributions in the region from ^ — 2ato^ + a, are equal to 3.5 and 
5.0 GeV, without and with overlaying events, respectively. 

The distributions of the Wcorr, obtained for the signal and background events, with over- 
laying events included, are shown in Fig. 7. The most precise measurement of the Higgs-boson 
cross section is obtained for the mass window Wcorr between 110 and 150 GeV, as indicated by 
arrows. In the selected Wcorr region one expects, after one year of the Photon Collider running 
at nominal luminosity, about 6900 reconstructed signal events and 8800 background events (i.e. 
S/B ^ 0.8). This corresponds to the statistical precision of: 



A 



T{h 77)Br(/i ^ bb) 



1.8%. 



T{h 77)Br(/i bb) 

We have performed a full simulation of signal and background events also for Mh = 130, 140, 
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Figure 6: Corrected invariant mass, Wcorr, distributions for selected 77 — /i — 66 events after 
6-tagging, for Mh = 120 GeV, obtained without and with overlaying events (OE). 

150 and 160 GeV choosing optimal e~e~ beam energies for each Higgs-boson mass. Statistical 
precision of T{h 77)Br(/i bb) measurement was estimated in each case. It is equal to 
1.8%, 2.1%, 3.0% and 7.1%, respectively. These results, together with described above result 
for Mh = 120 GeV, are presented in Fig. 8. For comparison our earlier results [8], obtained 
without overlaying events, are also shown. 



5 Conclusions 

We performed a realistic simulation of SM Higgs-boson production in the Photon Collider at 
TESLA, 77 — i> /i — i> bb, with the NLO background, corrections for escaping neutrinos and - for 
the first time - with the realistic 6-tagging and overlaying events. Our analysis shows that for 
Mh = 120-160 GeV the two-photon width of SM Higgs boson can be measured with a statistical 
precision 1.8-7.1%. 

The obtained accuracies are in a rough agreement with the results of a previous analysis, 
based on the idealistic Compton spectrum [12]. As shown in [6], the realistic photon-photon 
luminosity spectrum is more challenging for a precise determination of r{h — >■ 77)Br(/i bb). 
The precision of the measurement has been improved after applying a correction for escaping 
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Figure 7: As in Fig. 5, for the corrected invariant mass, Wcorr, distributions. 

neutrinos and a mass-window cut. In the present analysis we include overlaying events and, 
as expected, the measurement precision decreases. However, a realistic 6-tagging used here, 
resulting in a better 66-tagging efficiency than estimated earlier, allows to counterbalance this 
effect. 

The measurement discussed in this paper can be used to derive the partial width T{h 77) 
[1]. For example, for a Higgs boson of = 120 GeV we estimate a precision 1.8% on the 
measurement of T{h 77)Br(/i — bb). Assuming that Br{h bb) will be measured at the 
e+e~ Linear Collider with a precision 1.5% [15], the partial width r(/i — > 77) can be extracted 
with an accuracy of 2.3%. Using in addition the result from the e~^e~ Linear Collider for 
BT{h —>■ 77) [16], one can also extract Ftot with a precision of 10%. 

For higher masses of the SM Higgs boson other decay channels should be considered, see 
e.g. [17]. 



Acknowledgments 

We would like to thank A. De Roeck, K. Monig and T. Kuhl for valuable discussions. M.K. ac- 
knowledges partial support by the Polish Committee for Scientific Research, Grants 2 P03B 
05119 (2003), 5 P03B 12120 (2003), and by the European Community's Human Potential Pro- 
gramme under contract HPRN-CT-2000-00149 Physics at Colliders. 



10 



^ 8 



1-^ , 



T 



u 
.o 



• withOE 
T without OE 




115 120 125 130 135 140 145 150 155 160 165 



M, [GeV] 



Figure 8: Statistical precisions of r{h 77)Br(/i —>■ bb) measurements for the SM Higgs boson 
with mass 120-160 GeV with and without overlaying events (OE). The lines are drawn to guide 
the eye. 
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